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We add a Heisenberg interaction term oc A in the one-dimensional SU(2)®XY spin-orbital model 
introduced by B. Kumar. At A = 0 the spin and orbital degrees of freedom can be separated by a 
unitary transformation leading to an exact solution of the model. We show that a finite A > 0 leads to 
spontaneous dimerization of the system which in the thermodynamic limit becomes a smooth phase 
transition at A —> 0, whereas it remains discontinuous within the first order perturbation approach. 
We present the behavior of the entanglement entropy, energy gap and dimerization order parameter 
in the limit of A —> 0 confirming the critical behavior. Finally, we show the evidence of another 
phase transition in the Heisenberg limit, A — > 00 , and give a qualitative analytical explanation of 
the observed dimerized states both in the limit of small and large A. 
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I. INTRODUCTION 

A Peierls transition being a distortion of the periodic 
lattice of a one-dimensional (ID) crystal has been known 
for more than 80 years [1], Peierls’ theorem states that 
a ID equally spaced chain with one electron per ion is 
unstable to the dimerization, i.e., a more favorable con¬ 
figuration can be achieved if every ion moves closer to its 
one neighbor and further away from the other. This is 
a purely kinetic effect related with the nesting of the ID 
Brillouin zone. Such dimerization involves charge and 
lattice degrees of freedom but a similar process can be 
observed in the spin-orbital models [2, 3] where the in¬ 
teraction within the orbital sector dimerize the spin or¬ 
der. Such an effect was observed by a neutron scatter¬ 
ing in the perovskite vanadate YVO3 [4] as theoretically 
described by Horsch and coworkers [5]. Similarly, for 
an iridate triangular-lattice superconductor Ir 1 _ a ,Pt a; Te2 
the photoemission and model studies have revealed that 
an orbitally induced Peierls effect governs the charge and 
orbital instability in this compound [6]. A strict, ID spin- 
orbital analog of the Peierls’ dimerized state was found in 
NaTiS^Og [7], as shown by Raman scattering measure¬ 
ments, where the dimerized spin state at low tempera¬ 
ture is driven by the orbital fluctuations at high tem¬ 
perature. A quite different situation emerges when the 
orbital fluctuation do not drive the spins’ dimerization 
but the orbital order makes the electrons’ hopping quasi- 
ID so that the classical Peierls transition can occur [8]. 
Such an orbital-selective Peierls dimerization was found 
in the t 2 g spinel MgTi 2 04 by an optical measurements. 
On the other hand it was shown that the notion of the 
Peierls instability also applies to the strongly interact¬ 
ing electrons coupled to the phonons within the Holstein 


model [9] which allows to explain the orbital order in the 
undoped manganite LaMnOa [10]. 

Spin-orbital models are crucial from the point of view 
of frustrated magnetism [11-13] as the orbital interac¬ 
tions, being typically of lower symmetry than the SU(2)- 
symmetric spin interactions, can lead to frustration in 
both sectors even on a square lattice. This often hap¬ 
pens in the transition-metal oxides. When the 3d or¬ 
bitals are partly filled, electrons localize due to large 
on-site Coulomb interaction and superexchange between 
magnetic ions includes both spin and orbital degrees of 
freedom that are strongly interrelated [14]. The orbital 
degeneracy leads then in many cases to a dramatic in¬ 
crease of quantum fluctuations [15], which may trigger 
exotic order [16, 17], or may stabilize a spin-liquid [18, 19] 
when different states compete near a quantum critical 
point. In general these fluctuations make the ground- 
states of the spin-orbital models entangled, as it hap¬ 
pens in the archetypal Kugel-Khomskii model in one-, 
two-, three-dimensional and bilayer cases, respectively 
[16, 17, 20, 21). Similarly, in the 5=1/2 SU(2)0SU(2) 
chain [22, 23], both ground state [24] and excited states 
[25] are entangled and the 5=1 SU(2)<g>SU(2) chain 
which plays a prominent role in the vanadium perovskites 
[13, 26]. In exceptional cases can such ID models be 
solved exactly, for example at the SU(4) point [27] or for 
a valence-bond state [28] of alternating spin and orbital 
singlets similar to the Majumdar-Ghosh state in a ID J\- 
J 2 spin chain [29], but even in these situations the spins 
and orbitals cannot be separated from each other. 

In real materials the symmetry between spin and or¬ 
bital interactions is absent. Orbital interactions gener- 
ically have lower symmetry than spin ones [30], being 
usually Ising- or XY-like [31]. The XY case is quantum 
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and in general the orbitals cannot be separated from the 
spins [32]. Therefore the ID SU(2)(g)XY model intro¬ 
duced by Kumar [33] is very exceptional because; (i) it is 
exactly solvable and (it) by a change of basis, the 5 = 1/2 
spins decouple from the orbitals in an open chain. The 
orbital interactions remain formally unchanged but the 
spin ones are gauged away. The spins then appear free 
and the ground state has large degeneracy (2 L for chain 
length L) [33], 

Quite surprisingly, the properties of the Kumar model 
are determined by topology. When the chain is closed 
and the model becomes periodic the exact solution still 
exists, as shown in Ref. [34], but the spins can no longer 
be gauged away completely. Contrarily, the gauge is ac¬ 
cumulated on the closing orbital bond so that the whole 
spin sector enters the orbital problem only through the 
boundary condition, in analogy to Aharonov-Bohm mag¬ 
netic flux through the ring. This partially lifts the macro¬ 
scopic degeneracy of the ground state of the open chain 
and splits the manifold of the 2 L states into a multiplet 
of low-lying states. These are the states with topolog¬ 
ical excitations of a non-local nature having quadratic 
dispersion in the total quasimomentum carried by the 
spins [34]. Such a multiplet of states connected by the 
non-local excitations is a characteristic feature of another 
model originating from the orbital exchange, namely 
the two-dimensional orbital compass model (OCM) [35], 
which was considered, iter alia, in the context of topo¬ 
logically protected qubits [36]. From this point of view 
it was important to study how the small admixture of 
the SU(2) Heisenberg interactions, which is unavoid¬ 
able in the physical implementations employing arrays of 
Josephson junctions, will modify the multiplet structure 
in OCM [37], 

In this paper we will study the spin-orbital Kumar 
model with an extra antiferromagnetic (AF) Heisenberg 
term for the nearest-neighbor (NN) spins, namely the 
Kumar-Heisenberg (KH) model. As we will see further 
on, this model is no longer exactly solvable but the Ku¬ 
mar gauge transformation can be still applied. As a result 
the pure Kumar term gets completely or partially decou¬ 
pled depending on the topology of the system but the 
Heisenberg term is subject to the spin permutation that 
depends on the configuration of orbitals in their canoni¬ 
cal basis. In other words, under the transformation the 
NN spin interaction gets smeared out and this smearing 
is governed by the delocalization of the up/down orbitals 
caused by the quantumness of the XY orbital interac¬ 
tions. This situation is interesting because in this way 
the spin interactions can get frustrated - the smeared in¬ 
teraction remains AF but now the further neighbors are 
involved. On the other hand, as the pure Kumar model 
leaves the spins practically completely uncorrelated, we 
expect that any small admixture of spin interactions A 
will order the spins in some way. This is a favorable sit¬ 
uation as it will enable us to use perturbation expansion 
in the manifold of degenerate ground-states of pure Ku¬ 
mar model. To examine the region, where the coupling 


between the spins cannot be treated as a small pertur¬ 
bation, we use the density matrix renormalization group 
method (DMR.G) [38-42], which gives an accurate de¬ 
scription of the ground state in one dimension. This also 
provides an independent way to verify the perturbational 
results. 

The paper is organized as follows. In Sec. II we intro¬ 
duce the Kumar-Heisenberg (KH) model, in Sec. Ill we 
show the exact solution of the periodic KH model based 
on Ref. [34], in Sec. IV we apply this solution to the fer¬ 
romagnetic, A < 0, case and get the exact ground state, 
and in Sec. V we explain the spin dimerization within 
a first order expansion in a small A > 0 which is a valid 
description in the limit of short chains. In Sec. VI we 
explain why this limit cannot describe an infinite system 
where a second-order phase transition occurs, as shown 
in Sec. VII. In Sec. VIII we describe a dimerized state 
and give a mean-field mechanism of the orbital dimeriza¬ 
tion for large A. In this Section we also show the evidence 
of a second phase transition in the limit of A —> oo. The 
conclusions are presented in Sec. IX and in Appendix 
A we give further details on the effective spin couplings 
introduced in Sec. V. 


II. KUMAR-HEISENBERG MODEL 

The Hamiltonian of the KH model reads 

n = n 0 + v, (i) 


where 


L 

H 0 = jJ2 X M+1 (t+t- +1 +t-T+ +1 ) (2) 

1=1 

is the Kumar model and 

L 

V = \ y^yiai+i (3) 

1=1 


is the Heisenberg term. Here r’s are orbital Pauli matri¬ 
ces and Xi j = (1 +3i(jj) /2 is a spin transposition oper¬ 
ator, defined by another set of Pauli matrices a. In the 
case of periodic boundary conditions (PBC) o l+i = o\ 
and t l+i = T\. For open boundary conditions (OBC) 
<j l+ i = 0 and t l+1 = 0. 

The Hamiltonian can be transformed into base where 
T~Lq appears as a purely orbital contribution. For this we 
use the unitary operator 


l- 1 r 


"=II 


i —l 


1 - t, 


/-f-i 


+ 


1 + T, 


l~\~l , 




( 4 ) 


where Xi+i,i is a spin permutation operator composed of 
the spin transpositions AQj: 


Xl+i,i = Xt+ijXi'1-i...XS' 2 X 2 , 1 . 


( 5 ) 



3 


Now according to Ref. [34] the transformed Hamiltonian 
takes a form 

n'=u^nu=H' 0 + V, ( 6 ) 

where 

L- 1 

n 'o = J J2( T f R+i+ h - c ■)+j(r[ 1) Ri ) t^ rr+h.c.),(7) 
1 = 1 

and 

L 

V = \y^aiai+i. ( 8 ) 

l=i 

Here R ^ is a cyclic permutation of spins at sites l = 

( 2 ) 

1 ,...,TV by 1 site and R]~ is the same permutation of 
spins at sites l = TV + 1, L: 

R^ : {<ji, {<r 2 , • ••, <Uv, c?i}, 

R\ '■ {f?AT+l, - <Jl} —> {<Uv+ 2 , •••, f?L, &N+ 1 }- (9) 

Here TV is a good quantum number of ”up” (rf = +1) 
orbitals: 

n = Y, 1 -^- ( 10 ) 

1=1 

The primed spins S[ are entangled spin-orbital operators: 

<4 = |f 2 )a /( ^)(f z |, (11) 

with | t z ) = |rf, rf,... , r£) and a permutation of chain 
sites 

/a,**) = *+^(i+?f)(i-o+ E \( 1 + t p )• ( 12 ) 

p=;+i 

<j[ are the original spins d) after a permutation of the 
chain sites dictated by the orbital configuration in the 
rf-eigenbasis. 


III. SOLUTION OF THE KUMAR MODEL 

The Hamiltonian "Hq is exactly solvable [34]. On all 
bonds but the closing one it involves only orbital de¬ 
grees of freedom and even on the closing bond the spin 
part affects only the boundary conditions via the cyclic 
permutations R\ ’ . These can be diagonalized and 
turned to the phase factors e llCl ’ 2 with quasi-momenta 
(spin currents) K,\ = and Here 

n± = 0,..., TV — 1 and ?r 2 = 0,..., L — TV — 1. Then the 
diagonalization of "Hg can be completed by the Jordan- 
Wigner (JW) transformation of the form rf = 1 — 2ri/ and 
T; + = c; 11^(1-2^), where n; = qc/ and ci annihilates 


a JW fermion. The last step is a Fourier transformation, 
we have 

L—l 

Rq = J E c\ +1 ci+Je^ 2 ' ni ^c\c L + h.c. = 2 J^cj'c fc cosfc, 

1=1 k 

(13) 

where the phase 2i r$ = /Ci+/C 2 — 7 r(L — TV—1) is twist¬ 
ing the boundary condition, Cl+i = e 27 ri$ Ci, just like a 
magnetic flux $ through the periodic ring. The standard 
cosine dispersion of the fermions involves quasi-momenta 
k quantized as k = ^ (m + <!>) with m = 0,..., L — l. 

As shown in Ref. [34], for PBC the ground state (GS) 
of the model is the Fermi sea of the JW fermions at half¬ 
filling, 

TV = L/2, (14) 

with a zero spin current, 

/Cr+/C 2 =0. (15) 

Thus the GS wave function of Hq factorizes into a prod¬ 
uct of an orbital and spin state: 

i°) = 4 , • • • 4 W \ vac )c ® i^i) (1) i-^i) (2) , (i 6 ) 

with |uac) c being a vacuum of the JW fermions. The 
orbital state is a Fermi sea with occupied quasimomenta 
kj = -j^(—j + g) — f for j = 1,..., TV such that cos kj < 0. 

The spin states |/C ) < + 2 ' 1 are the /C-momentum eigenstates 
of the operators with an eigenvalue e l,c . The 

ground state is degenerate with respect to the TV values 
of/Ci. 

For OBC the spin state is completely arbitrary and the 
orbital state is a half-filled Fermi sea with OBC. Since in 
this paper we are interested in the thermodynamic limit, 
from now on we focus on OBC. 

The Kumar model is the special case of the Heisen¬ 
berg coupling A = 0 in Eq. (3). Its 2 L -fold degenerate 
ground state is very sensitive to the perturbation. In the 
following we scan different regimes of A. 


IV. FERROMAGNETIC COUPLING: A < 0 

This case is easy to solve in the Kumar basis, where 
the ferromagnetic spin state is the ground state of both 
the Heisenberg term (8) and the Kumar model (7). The 
orbital ground state is a half-filled Fermi sea. 

Back in the physical representation, permuted spins 
are still ferromagnetic, the exchange operators X{j = 
1 , and the orbitals are in the ground state of the anti¬ 
ferromagnetic XX-chain, see Eq. (2). The ground state 
is a product of the ferromagnetic spin state (FM) and 
the orbital state - there is no spin-orbital entanglement. 
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orbitals up 



Figure 1. The schematic view of the bonds connecting the 
two spin subsystems (upper and bottom legs of the ladder) 
in the first-order Hamiltonian in Eq. (17). The spin bonds 
connecting site l = 4 of the orbital-up subsystem with sites of 
the orbital-down subsystem are marked. Their color intensity 
increases with the strength of the coupling J ij. 


V. PERTURBATIVELY WEAK COUPLING: 

0 < A < 1/L 

On the antiferromagnetic side the solution is no longer 
so simple. In order to get a rough idea how the highly 
degenerate Kumar ground state may respond to a weak 
Heisenberg perturbation, it is tempting to try a pertur¬ 
bation theory in weak A. In the Kumar representation for 
OBC, the spin state is arbitrary and the orbital state is 
a half-filled Fermi sea with a gap ~ 1/L. When the per¬ 
turbation is weaker than the orbital gap, 0 < A <C 1/L, 
one can treat the spins with the degenerate perturbation 
theory. 

To first order in A, the orbital state does not change, 
protected by the orbital gap, and the spin degeneracy is 
removed by a Hamiltonian 

L 1 

K = <V'> = A Y, 2 J ^ % ( 17 ) 

i,i =i 

where the average (..) is taken in the orbital Fermi sea. In 
the physical representation only nearest-neighbor (NN) 
spins are coupled. 

After the Kumar transformation, each physical NN- 
coupling contributes to (is smeared over) many .l hl (see 
Appendix A for more details). The transformation is a 
permutation that maps the N consecutive spins on empty 
sites (orbitals up) to spins i = N,.... 1 and the N con¬ 
secutive spins on occupied sites (orbitals down) to spins 
i = N + 1, ...,21V. The exact matrix J 2 J is described in 
more detail in Appendix A. It turns out to be dominated 
by the couplings between the two IV-spin subsystems, as 
shown in Fig. 1. The strongest bonds are along the 
rungs of the spin ladder. When we keep only them, then 
the ground state will become a product of spin singlets 
along the rungs. We compared the energy of this trial 
state to the energy from TL\ with the full matrix J. For 
L = 16 the energies differ by less than 0.00004 in case of 
A = 0.01 and by slightly more than 0.0003 for A = 0.1. 

A full trial state is the spin trial state times the or¬ 
bital Fermi sea. It can be transformed back to the phys¬ 
ical representation. For L = 16 the result shows strong 
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Figure 2. Finite-size scaling results for the block entropy of 
an open chain as a function of 1/L (up to length L = 600) for 
several values of A. For each A the entropy saturates for small 
enough 1/L where the thermodynamic limit is achieved. 

dimerization of spins i.e., (<7;<7/+i) < 0 on odd NN bonds 
and = 0 on even ones. The odd bonds take large val¬ 
ues: {-2.590, -2.183, -2.063, -2.018, -2.018, -2.063, 
—2.183, —2.590}. This perturbative result uncovers a 
tendency of the system towards spin dimerization. The 
transition to the dimerized state is discontinuous (first- 
order) by the very nature of the degenerate perturbation 
theory. 

However, the degenerate perturbative treatment is not 
justified when A becomes comparable or stronger than 
the orbital gap ~ 1/L, hence it must fail, at least quanti¬ 
tatively, in the thermodynamic limit. This is why in the 
following we use DMRG for OBC to obtain numerically 
exact results in this limit. 


VI. CROSSOVER TO THE THERMODYNAMIC 
LIMIT: A ~ 1/L 

In our DMRG calculation we used the dynamic block- 
state selection algorithm. [43, 44] The a priori value of 
the quantum information loss was set to x = 10 _ 7 , which 
required to keep block states up to 4000. Typical trun¬ 
cation errors were in order of 10 -8 . We considered finite 
chains with OBC up to length L = 800. 

The thermodynamic limit is reached when the system 
size is longer than its correlation length £. The crossover 
to this regime can be seen e.g. in the entanglement (von 
Neumann) entropy, s(A, L/2), between the two halves of 
the chain. [45-47] This quantity can be obtained from the 
appropriate reduced density matrix, p L / 2 (A): 

s(A, L/2) = -Trp i/ 2 (A) lnp i/ 2 (A), (18) 
where Pl/ 2 (A) denotes the reduced density matrix of the 
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Figure 3. The entanglement entropy of a half-chain as a func¬ 
tion of A in the thermodynamic limit. Note that the rr-scale 
is logarithmic. The points are the DMR.G results, while the 
solid line is the best fit using Eq. (20). 

subsystem containing sites 1, ...,L/ 2. In Fig. 2, we show 
the entropy for different system sizes and small A <C 1. 

For small system sizes, or large 1/L, the entropies 
for different A collapse. This is the regime of validity 
of the degenerate perturbation theory, where the actual 
strength of a weak A <C 1/L is irrelevant. 

For each A the entropy is growing with the system size 
L until it saturates below a crossover value of 1/L. This 
saturation is characteristic for the thermodynamic limit. 
Interestingly, the weaker A is the longer L is needed for 
the saturation, hence the correlation length £ must in¬ 
crease or even diverge when A —> 0. In the following Sec¬ 
tion, we consider mainly ‘saturated’ results in the ther¬ 
modynamic limit for all physical observables. 

VII. THERMODYNAMIC LIMIT: 
SECOND-ORDER TRANSITION FOR A —> 0+ 

Firstly, we consider the behavior of the block entropy 
of a half-chain, s(A), in the thermodynamic limit: 

s(A) = lim s(A, L/2). (19) 

L—>o o 

This is shown in Fig. 3, where x-scale is logarithmic. 

The entropy diverges as we approach the critical point. 
Near the critical point at A —► 0 + the entropy is expected 
to behave as [47] 

s(A) = ^ln(O = -yln(A/A 0 ), (20) 

where c is the central charge and £ = (A/Ao) - " the cor¬ 
relation length with the exponent v. The best fit is 
A 0 = 3.9 ± 0.5 and cv = 1.11 ± 0.02. 



1/L 


Figure 4. The finite-size scaling of the spin gap for several 
values of A. The solid lines denote the result of our fit using 
Eq. (21). 

The diverging block entropy indicates a diverging cor¬ 
relation length / in the system and a vanishing gap 
A = A 0 \ zv , where z is the dynamical exponent. There¬ 
fore we investigated the first few low-lying excitations in 
the model. We found that the first excitation for A <C 1 
is a spin flip, that is, the smallest gap is given by the 
spin gap. To obtain the bulk values of the gap, A(A), 
we performed extrapolation to the thermodynamic limit 
with the fitting ansatz 

A(A, L) = A(A) + A/L + B/L 2 , (21) 

where A(A), A and B are free parameters. This ansatz is 
motivated by the fact that the finite size corrections are 
expected to be algebraic for OBC. The finite-size scaling 
procedure is demonstrated in Fig. 4. 

The extrapolated gap as a function of A is shown in 
Fig. 5 with a log-log scale. 

The data points can be fitted well with a linear, that 
is, the gap exhibits a power-law dependence on A: 

A(A) = A 0 A W , ( 22 ) 

where Ao = 3.12 ± 0.02 and zv — 1.088 ± 0.002. This 
means that the gap vanishes continuously in contrast to 
the naive prediction of the perturbation calculation, and 
a second-order phase transition occurs at A = 0 + . 

The same applies for the dimerization, D, which is 
defined as 

D(A) = lim |s(A,L/2) — s(A,L /2 — 1 )|. (23) 

L—> oo 

This is shown in Fig. 6 with a log-log scale. The dimer¬ 
ization - similarly to the gap - behaves as 

D( A) = D 0 \P, (24) 

where D 0 = 1.22 ± 0.03 and /3 = 0.204 ± 0.004. 
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A 

Figure 5. The spin gap extrapolated to the thermodynamic 
limit as a function of A. Note that both scales are logarithmic. 
The points are the DMRG results, while the solid line is the 
best fit using Eq. (22). 




Figure 6. The extrapolated dimerization order parameter, D, 
as a function of A. The points are the DMRG results, while 
the solid line denote the fit using Eq. (24). Note that both 
scales are logarithmic. 


VIII. DIMERIZATION FOR FINITE A BOTH 
LARGE AND SMALL 

The previous Section focused on the critical behaviour 
in the limit A —> 0 + . This is a second-order phase transi¬ 
tion, where the symmetry broken by dimerization is re¬ 
stored. In this Section we paint a general picture of the 
dimerized phase at A > 0. We distinguish two regimes, of 
large and small A, where the nature of the dimerization is 
qualitatively different. They are separated by a smooth 
crossover. 

It is interesting to analyze the behavior of the block en- 
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A 

Figure 7. The extrapolated block entropy, s(A), and dimer¬ 
ization order parameter, D( A), as a function of A are shown 
in panels (a) and (b), respectively. 


tropy and the dimerization order parameter for larger A 
before going into the details of the dimerized phase. They 
are shown in Fig. 7. It is seen that around A ~ 0.7 the 
block entropy and the dimerization order parameter ex¬ 
hibits a local minimum and maximum, respectively. This 
can be understood as follows: when the dimerization is 
the strongest the ground state is the less entangled, which 
manifests in the minimum value of the block entropy. We 
address what happens for A —>■ oo later in this section. 


A. Numerical findings 

Numerical observations are as follows. For any A we 
observe both spin and orbital dimerization in the sense 
that the NN correlations alter between large/small values 
for odd/even bonds as shown in Fig. 8. The difference be¬ 
tween spin and orbital sector is that the spin correlations 
remain always AF while the orbital ones can alternate 
between positive/negative values for small A and remain 
FM only for large A. The behavior of the long-range 
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Figure 8. DMRG results for the NN spin (panel (a)) and 
orbital (panel (b)) correlations for the open chain of L = 140 
sites as a function of the site number l for A = 0.01 and A = 5. 


spin and orbital correlations is different in the small and 
large A regime. While the correlation (aiai+n) always 
decays exponentially with R it remains constantly neg¬ 
ative for small A while for large A it alternates between 
negative/positive values. Thus we say that for small A 
we have strong dimerization of the spins. Because of 
the negative sign of the correlations this is a spin-singlet 
dimerization. The behavior of the orbital correlations 
seems to be complementary, the long range correlations 
(t ; + t ; ^ r ) are constantly positive for large A, while for 
small A they alternate with R. In both regimes they de¬ 
cay exponentially with R. Thus, analogically to the spin 
sector, we say that for large A the orbitals exhibit strong 
dimerization in the orbital-triplet state. Summarizing, 
we have found strong spin-singlet dimerization for small 
A and strong orbital-triplet dimerization for large A. 

Finally, we note that the transition between the two 
regimes, of strong dimerization of either spins or orbitals, 
is a crossover without any features of a phase transition. 
However, one can easily observe a change in the behavior 



Figure 9. DMRG results for the spin and orbital correlations 
for the open chain of L = 140 sites as a function of the site 
number l for A = 0.01 and A = 5. The vertical axis has a 
logarithmic scale. 
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Figure 10. Finite-size scaling results for the block entropy of 
an open chain as a function of 1 /L (up to L = 800) for several 
decreasing values of J and A = 1. For each J the entropy sat¬ 
urates for small enough 1/L where the thermodynamic limit 
is achieved. 


of the entanglement entropy of the half-system s(A) in 
Fig. 7 (a); while it is rapidly decreasing for small A, it 
becomes slowly increasing for larger values of A with a 
shallow minimum at A ~ 0.7, where the dimerization is 
the strongest. 

At the same time, in addition to the qualitative change 
of the spin correlation triggered by increasing A, there is 
a big difference in their correlation lengths. They are 
shown in Table I. From the semilogarithmic plot of Fig. 
9 we can learn that for small A both correlations decay 
with almost the same exponent, while for large A the spin 






































































correlation function has much larger correlation length 
than that of the orbital correlation function. This sug¬ 
gests that there is a second phase transition in the system 
located at A —> oo or J = 0. This is not totally unex¬ 
pected as in this limit the KH model becomes simply the 
Heisenberg model in the spin sector which is known to 
be critical. This observation is confirmed by the behavior 
of the block entropy of the half system, s(L/ 2, J), which 
exhibits similar scaling properties as J —> 0 for A = 1, see 
Fig. 10, to what we have seen for A —► 0 in Eq. (20). In 
the present case we fit the product cv = 1.85±0.01. Since 
c = 1 for the Heisenberg chain, we obtain an estimate for 
the correlation length exponent v = 1.85 ± 0.01. 


A 

^orbital 

£spin 

0.01 

5.52 

4.65 

5 

4.72 

18.18 


Table I. The correlation lengths of the spin and orbital corre¬ 
lation functions for several values of A. 


B. Qualitative explanation 

The dimerization of spins involving constantly nega¬ 
tive long-range correlation function can be understood 
in terms of the first order perturbative expansion of Sec. 
V and especially looking at Fig. 1. As we can see, in the 
Kumar basis the singlet bonds are mainly on the ladder’s 
rungs connecting spins at sites N + i and N — i + 1. In 
the physical basis the rung singlets become, in first ap¬ 
proximation, spin singlets on every second bond. In this 
approximation the correlation function is strongly nega¬ 
tive only within a given singlet and zero anywhere else. 
However, when we take into account fluctuating positions 
of the up- and down-orbitals in the orbital Fermi sea, 
then the singlets on every second bond become smeared 
over neighboring sites. This smearing implies that: (z) 
singlets placed every second bond are not perfect, i.e. 
0 > ((j/cq+i) > —3 and (zz) there are negative spin cor¬ 
relations between sites l and m £ [l — 61, l + 61] where 
61 is the legth scale at which a single up orbital is de¬ 
localized within the Fermi sea. In this framework the 
strongly spin-dimerized ground state of the KH model 
can be understood as a spin liquid of strongly resonating 
singlets. 

The triplet dimerization of the orbitals for large A can 
be easily understood in the mean-field way similarly as it 
was done in Ref. [2], In this limit the leading part of the 
Hamiltonian is a pure Heisenberg term V for the spins. 
Thus we can use a mean-field decoupling of a spin-orbital 
term in 7-l 0 : 

x i,i+i( T f T r+i+ T r t i+ i) ~ ( x u+i) (r + r+ i+Tfr+i) + 

X u+1 <R + r / ” 1 +T,"r+ 1 ) - ( x i,i+i) (r+T^+Tf 7-,+ i) . 

( 25 ) 


Kumar 


FM 


i 


dimerized 


A = 0 


Heisenberg 

-► A 

A = oo 


Figure 11. Schematic phase diagram of the Kumar- 
Heisenberg model for J = 1. 


The spins are AF with (a/tT;+i) < —1, hence for each 
bond we have (Xu + 1 ) < 0 and, according to the first line 
of Eq. (25), the orbitals order ferromagnetically. The in¬ 
stability towards dimerization, meaning here formation 
of the orbital triplets, can be demonstrated by imagining 
a self-consistecy loop for a non-uniform MF Hamiltonian. 
If, for any reason, a bond has increased the orbital cor¬ 
relation, then, according to the second line of Eq. (25), 
the spins will feel an increased tendency towards anti¬ 
ferromagnetism. Then again, the first line of Eq. (25) 
implies that the orbital bonds will be increased and so 
the spin bond will follow. Such a loop clearly leads to 
maximalization of the orbital correlations at some bonds 
and, since we cannot have a triplet at every bond, we 
must have the dimerization of the orbitals. We verified 
this argument by solving numerically a non-uniform MF 
Hamiltonian for L = 16. 


IX. CONCLUSION 

We have found that an admixture of the Heisenberg 
interactions A between the spins causes a spontaneous 
dimerization of the Kumar model. The dimerization is 
both in the spin and orbital sector but it becomes strong 
for the spins in the regime of small A and for the orbitals 
in the opposite limit. We have shown that the dimeriza¬ 
tion of spins can be understood by means of the Kumar 
transformation, i.e., in a basis where the spins and or¬ 
bitals in the Kumar Hamiltonian get decoupled. On the 
other hand we show that the dimerization of orbitals for 
large A can be understood be a mean field mechanism in 
the physical basis. 

We have seen that the perturbative expansion in the 
Kumar basis, however useful for understanding the spin 
dimerization, fails to capture the second-order phase 
transition found by DMR.G in the thermodynamic limit. 
Whereas the perturbative prediction is always a discon¬ 
tinuous transition, the numerical study of the system of 
sizes up to L = 600 allowed us to observe the divergence 
of the block entropy, vanishing of the gap and the dimer¬ 
ization order parameter in the regime where we approach 
the transition point at A = 0 slower than the gap in the 
pure Kumar model closes, i.e., in the regime where the 
perturbative approach must fail. Finally, by observing 
that the spin correlation length is strongly increased as A 
grows we have deduced that there is another phase tran- 
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sition in the limit of A —» oo (or equivalently J = 0). This 
we have confirmed by the finite-size scaling of the block 
entropy for decreasing J. The phase diagram is shown 
schematically in Fig. 11. 

Finally we note that however the full dimerization in 
orbital space can be easily understood by the mean- 
field approach, similarly as it was done for YVO 3 [2] 
at finite temperature, but here in zero temperature, 
the spin dimerization is more complex and its mecha¬ 
nism can be revealed only by the Kumar transforma¬ 
tion. This is a novel feature that was not found in the 
similar spin-orbital models, namely SU(2)(g>XY [31] and 
SU(2)®SU(2) [31, 48] one, where in principle the dimer¬ 
ization could be captured by a simple variational wave 
functions in the physical basis and confirmed via DMRG. 
We argue that such a novel dimerization can be observed 
in the (quasi) one-dimensional compounds with active 
spin and orbital degrees of freedom. 
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Appendix A: The effective spin couplings 

The effective, first-order couplings between the spins 
can be expressed as, 

N 

= P(p) (Al) 

Pl,---,PN Z = 1 

where the first sum is an average over a probability dis¬ 
tribution 

r®= n (A2) 

for empty sites {pi, ...,pjv} in a half-filled Fermi sea. Here 
we consider periodic boundary conditions to minimize 
boundary effects. For N = L/2 the permutation maps 
the N consecutive spins on empty sites (orbitals up) to 
spins N ,..., 1 and the N consecutive spins on occupied 
sites (orbitals down) to spins N + 1,..., 2 N. In the phys¬ 
ical representation only nearest-neighbor spins are cou- 


( a ) 4 8 12 16 




Figure 12. Plots of the J coupling matrices of the effective 
spin Hamiltonians Hi for system sizes; (a)— L = 16 and 
(b)— L = 100. Colors tending to red mean high positive 
values, white fields mean zero values. Dashed grid is a guide 
for the eye. 

pled. After the transformation U each physical NN cou¬ 
pling (PNNC) contributes to (is smeared over) many J ij. 

The matrix J,;j is shown in Figs. 12(a) and 12(b) as 
color plots for a small and large L respectively. The ma¬ 
trix has two diagonal and two off-diagonal N x N blocks. 
The diagonal block i,j = 1,..., N originates from PNNC’s 
between spins on empty sites. Since only a pair of con¬ 
secutive empty sites can happen to be NN’s, the block 
couples only NN’s. The weakest coupling Jq^v originates 
from the relatively rare situation when the first and last 
empty sites happen to be also respectively the first and 
last sites of the chain, i.e., they are not only NN’s (across 
the closing bond) but also each of them occupies a defi¬ 
nite position. Jqjv = 0 for an open chain. The diagonal 
block i,j = N + 1,2 N originates from PNNC’s be¬ 
tween spins on occupied sites. In our half-filled orbital 
Fermi sea the two diagonal blocks are the same. 

The two identical off-diagonal blocks originate from 
PNNC’s between pairs of empty and occupied sites. The 
top values of the coupling are between the first and the 
last site: Ji,l ~ 0.705 for L = 16 and J 1 .l — 0.783 
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for L = 100. This strong coupling originates from the 
PNNC between the IV-th empty and the IV-th occupied 
site that are very likely to be next to each other. The 
coupling Jpjv+i originates from the PNNC between the 
IV-th empty and the 1-st occupied site that is possible 
mainly across the closing bond. This coupling is much 
weaker for an open chain where the IV-th empty and 1-st 
occupied sites are rather unlikely to be NN. However, the 
dominant feature of each off-diagonal block is its antidi¬ 


agonal belt. A coupling J N—i,j within this belt originates 
from a PNNC between i-th empty and j-tli occupied site 
with i close to j. These sites are most likely to be NN 
when their difference \i — j\ is not much greater than a 
variation of a position of the z-th empty (or the j’-th oc¬ 
cupied) site in the Fermi sea. This variation limits the 
width of the antidiagonal belt. 

We found that truncating the full matrix .L (jJ to its 
dominant antidiagonal belt is a very good approximation. 
We are assuming this approximation in the main text. 
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